Transient polarization dynamics in a CO2 laser 
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Abstract 

We study experimentally and theoretically the polarization alternation during the switch-on 
transient of a quasi-isotropic CO2 laser emitting on the fundamental mode. The observed transient 
dynamics is well reproduced by means of a model which provides a quantitative discrimination 
between the intrinsic asymmetry due to the kinetic coupling of molecules with different angular 
momenta, and the extrinsic anisotropies, due to a tilted intracavity window. Furthermore, the 
experiment provides a numerical assignment for the decay rate of the coherence term for a CO2 
laser. 
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I. INTRODUCTION 



Laser dynamics is commonly studied considering the electric field as a scalar variable, 
since in most systems the polarization state is imposed by anisotropies of the cavity. For 
instance, Brewster windows or gratings, generally used in gas lasers to close the laser tube or 
to select a vibro-rotational transition, impose a linearly polarized state of the laser emission. 
However, in perfectly cylindrical laser cavities without any elements to select a preferred 
polarization, the study of the dynamics includes the necessity of considering the vector 
nature of the electric field. 

Several theoretical works have been devoted to the study of the polarization dynamics 
of the quasi-isotropic laser, showing the important role played by the material variables. In 
particular, the degeneracy of the angular momentum states of the laser transition sublevels 
has been considered as the coupling source between different polarization states. Initial 
studies considered only stationary solutions |], but more recently dynamical models have 
been developed to explore the role of anisotropy due to the laser medium, that from now on 
we will call intrinsic anisotropy [|3], |], 0, ||. These models predict a rich dynamics even for 
the simplest transition, J = 1 — > J = 0. Later versions of these models include also linear 
and circular cavity anisotropies that we call extrinsic || |9| , or different level structure [JTO] . 



However, up to now few experiments have been performed on this subject. Experiments 
carried out on noble gas lasers reveal that in some systems it is necessary to consider the 
dynamics of the matter variables to fully understand the polarization features of an isotropic 
system |2], |ITTj. In other cases, the observations could be explained by a nonlinear coupling 



of the modes and residual cavity anisotropies [13, 14 



In this work we present experimental and theoretical results for the polarization com- 
petition dynamics in the transient state of a quasi-isotropic low-pressure CO2 laser. We 
show that the polarization dynamics during the switch-on is well described by means of a 
model including optical coherences (intrinsic anisotropy) and extrinsic linear anisotropies. 
These results extend our previous work in which the coherences were simplified 
metric cross coupling between the matter polarization and electric fields jnj , and provide a 
quantitative comparison between the amount of anisotropy due to the coupling of molecular 
angular momenta and that induced by an intracavity window. Furthermore, the decay rate 
of the coherences, defined in Ref. H], is shown explicitly to be closer to the decay rate of 
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the population inversion rather than to the decay rate of the induced polarization. 



II. EXPERIMENTAL RESULTS 

The experiment has been performed using an unpolarized Fabry-Perot cavity as shown 
in Fig. |l|. A total reflective flat mirror (Mi, reflectivity Ri=l) and an outcoupler mirror 
(M 2 ) with a reflectivity i? 2 = 0.914 set the cavity length at L =75 cm. A piezo translator 
(Pzt) is used to select the P(20) laser emission line and to adjust the laser detuning. 

The active medium, a mixture of He (82%), N 2 (13.5%) and C0 2 (4.5%) at a pressure 
of 25 mbar, is pumped by a DC discharge fixed at 6.1 mA when the threshold current is 3 
rnA. 

As we are interested in the transient dynamics, an intracavity chopper is used to induce 
a switch-on event at a repetition rate of 200 Hz. 

In order to control the linear anisotropies of the cavity, a ZnSe anti reflection coated 
window (W) has been introduced, which can be tilted accurately as shown in Fig.|I]. 

The polarization state of the laser emission has been analyzed by means of a wire grid 
polarizer (Pol), which has the property of reflecting one linear polarization of the incident 
radiation and transmitting the orthogonal one, with an extinction ratio of 1:180. The 
reflected and transmitted parts of the beam are directed to two HgCdTe fast detectors 
(Di, D 2 , 100 MHz bandwidth), whose sensitive areas (10 4 fim 2 ) are much smaller than the 
beam size. Both signals are recorded on a digital oscilloscope (Lecroy LT423L) with 500 
MHz bandwidth. 

In the stationary regime, we observe that the laser has two possible linear polarization di- 
rections which are orthogonal as far as we can measure. In Fig. |1] these cavity eigendirections 
are called H and V respectively. 

The window W can rotate by an angle 9 around the axis V. When 8 = 0, the polarization 



direction is determined by the detuning of the cavity, as usual 0, [| jlB] . This is illustrated 
in Fig.[2](a), where we plot the laser intensity along both eigendirections when the detuning 
is varied around the line center. A residual hysteresis around the transition is masked by the 
limited resolution of Pzt. At the transition point near the line center there is not a preferred 
polarization direction and the laser flips spontaneously from one to other as shown in Fig. 
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Our goal is to study the polarization behavior in this bistable region. In particular, we 
analyze of the switch-on transient state, in which the total intensity presents relaxation 
oscillations. In our system the pump strength imposes a relaxation frequency of 55 KHz. 

When the two polarized components are separated, they show oscillations in relative 
antiphase, which do not appear on the total intensity (Fig.^). The amplitude of these 
oscillations depends on the angle between the polarizer axes with respect to H-V, reaching 
a maximum when the analysis is performed at 45°. 

In Fig.|](a) we report an example of these polarization oscillations, when the cavity is 
tuned at the center of the line. The total laser output intensity (thick solid line) is displayed 
together with its two orthogonally polarized components analyzed at 45° (thin solid and 
dashed lines). 

These oscillations are always damped until they disappear, with a damping rate depend- 
ing on the cavity detuning. Precisely the oscillations are more persistent the closer is the 
detuning to the bistable region. This fact points to the competitive origin of the oscilla- 
tions. In Fig.[| (c) we show an example of the transient dynamics when the cavity detuning 
is slightly moved from resonance. It is worth to note that the system is so sensible to noise 
that the duration of the oscillations suffers slight variations, without apparent changes in 
the conditions. 

Successively, we studied the response of the system to small linear anisotropics driven by 
the tilt angle 9 of the intracavity window W around the vertical axis. 

Several parameters of the system are simultaneously affected by this action. As a conse- 
quence of the common increasing of the optical path for both polarization eigendirections, 
the line center moves with respect to the isotropic condition 9 = 0. In order to keep the 
bistable condition, the cavity length is varied to recover the resonance. 

In Fig.|2|((b) and (c)) we show the intensity of both polarization components for different 
9 values, showing that the gain profile is not significantly distorted with respect to the 
isotropic condition 9 = (Fig.0(a)). The detuning anisotropics are measured to be smaller 
than 5% of the total change. 

The other effect of the tilt angle 9 is to increase the cavity losses. In Fig.^ (b) we report 
the dependence of the laser intensity on the tilt angle 9 which is in agreement with the 
expected theoretical behavior J/J. Moreover, since the gain profiles for the two polarization 
directions do not change due to the tilt angle 9 (Fig. |2|), we can assume that the induced 
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loss anisotropics remain a small effect ( ~ 1%) compared to the total loss change. The effect 
of the losses anisotropy can be observed in the residual modulation in the total intensity 
(Figs.| (a),(c))- 

The dynamical effect observed when 6 is increased is the growing of the frequency of the 
polarization oscillations. Measurements made near the resonance condition show that the 
frequency rises as the tilt angle increased(Fig.|5] (a)). It can be seen that a tilt angle of 12° 
is enough to rise the frequency from 100 KHz to 400 KHz. 



III. MODEL AND NUMERICAL RESULTS 



Our theoretical approach is based on the theory of the isotropic laser developed in Ref. 
H where the optical coherences between upper levels are considered with the. This theory 
was developed for the simplest case (J=l — > J=0), while the transition involved in our 
system is much more complicated (J=19 — > J=20). However, this theory has been showed 



to predict also the behavior of lasers with a different level structure, as shown in Ref. [jTT 
Only first order coherences (Am = ±1) will be considered. Therefore, independently of the 
number of sublevels, there are only two kind of possible transitions, which generate a split 
of the population in two ensembles, in such a way that an anisotropy is induced in the active 
medium §. Furthermore we introduce an extrinsic linear anisotropy as done in Ref.||. 

The field will be decomposed in a circularly polarized base. Just losses and linear de- 
tuning anisotropies along the principal axes of the system will be included, but not circular 
asymmetries since our system does not show signs of dicroism. It reads as 

E R = k(P r - E R ) + iSE R - (a + i/3)E L , 
E L = k{P l - E L ) + i5E L - (a + if3)E R , 
Pr = -i±[Pr-D r Er-E l C], 

Pl = -1±[Pl ~ D L E L - E R C% (1) 
C = - 1c C-^(EIP r + ErP* l ), 

Dr = -7|| [D R -r + \{E R P* R + E R P R ) + \{E L Pl + E* L P L % 

Dl = -1\\[D L - r + \{E L Pl + E* L P L ) + \{E R P R + E R P R )) 

where ER(r,t), ^(f, t) are the slowly varying electric fields in the circular basis. 



Pn(r,t), Pn(r,t) stand for the matter polarization fields, Dn(r,t) and Di{f,t) are the re- 
spective population inversions. The field C(r,t) represents the coherence between the upper 
sublevels. We recall that, in a density matrix treatment, the polarization correspond to off 
diagonal matrix elements between upper and lower level of the radiative transition, whereas 
C is proportional to the of diagonal matrix elements coupling different angular momentum 
states of the upper level ]TT|] . The parameter 5 represents the detuning between the cavity 
and the atomic transition frequencies. The parameters a = (kv—Kh)/2 and (3 = — 5y)/2 
represent respectively the linear anisotropics in the losses and detuning with respect to the 
cavity H-V axes, where Ky, k>h, are the losses in the horizontal and vertical axis, and $y, 
5h are the corresponding detunings. 

In our low pressure CO2 laser, the polarization decay is 7^ = 4.4 ■ 10 8 s _1 and the 
inversion decay rate as 7y = 1.95 • 10 5 s _1 [[HJ. In view of numerical solutions, there is no 
reason to keep dynamical equations for variables whose decay rates differ by several orders 
of magnitude. Depending on whether 7 C ~ 7y or j c >> 7y, we can respectively reduce 
the number of differential equations from 7 to 5 (for E#, E/,,C, and D^) or 4 (for Er, 
El, D# and D^). In the first case, Pr, Pl are eliminated adiabatically and their values 
correspond to the stationary solutions of the third and fourth of Eqs. [I]. In the second case 
(fast C), also the C variable is eliminated adiabatically and it follows the relatively slow 
variations of the field and the population inversion. Whether to choose between five or four 
equations depends on the numerical value of j c . The 7 C parameter represents the coherence 
decay rate, whose value should be chosen between 7j_ and j\\ [f|. However, this parameter 
can not be directly measured, and it will be used as control parameter in order to fit the 
theory to the experimental results ||. We find that the optimal value is 7 C ~ 7y in all cases, 
which is also consistent with the observation that just linearly polarized states are found 
in the experiment. Indeed, a higher value of 7 C would give rise to a periodic modulation 
of the total intensity |J which has never been observed in our experiments. Hence, the 
coherence C decays on the same time scale as the population inversion and we must keep 
the corresponding time dependent equation. Thus our experimental system is dynamically 
modeled as a class B laser ||12j| , and Eqs. (1) reduce to the following system of two algebraic 
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equations: 



P R = D R E R + E L C 

P L = D L E L + E R C* (2) 



and five differential equations: 



E R = 


k(D r - 1)E R + i5E R + (kC - 


-(a + iP))E L , 




E L = 




- \)E r + i5E l + (kC* 


-(a + i/3))E R , 




C = 


-TiKca+^r+i^i 2 )) 


-±E R Ei(D L + D R )), 




D R = 


-7|| 


D R {l + \E R \ 2 )-r + ^ 


^D L \E L \ 2 + ^(E R E* L C*- 


H E* R E L C)) 


Dl = 


-7|| 


D L {\ + \E L \*)-r + \\ 


f 2 3 
Dr\E r \ + -(E R E* L C* - 





(3) 



Eqs. (0) show explicitly that the cross field coupling is due partly to a molecular coupling 
J = 1 — >• J = (intrinsic) and partly to the intracavity window W (extrinsic). Fig. |5| 
provides evidence that the main responsible for the alternation is the intrinsic anisotropy 
rather than the extrinsic one. Figure |5| (a) refers to the case where both coherences and 
linear anisotropics are acting. In Fig. [| (b), the extrinsic anisotropy has been removed, 
now the oscillation competition is not damped, but the alternation of the two polarization is 
maintained. Finally, in Fig. |5] (c) just the extrinsic anisotropies are acting; the competition 
does not appear any more. 

From the experimental observations we know that the most important effect of the window 
tilting is the increase of the total loss, while the losses and detuning anisotropies represent a 
minor effect. In addition, as the excitation current does not change during the experiment, 
the increasing of the losses affects the gain to loss ratio. As a consequence, the total losses 
K and pump parameter r are functions of 9 as follows 0: 

r(0) =r (l -a r sin(fl) 4 ) (4) 
K (6) = K (l + a k sm(9f) (5) 

where r Q = 2 and k = —^log(R) = 4.1 ■ 10 6 s _1 are the pump strength and total losses 
when the tilt angle 6 = 0. Here c is the speed of light, L = 0.75 m the cavity length and 
R = y/RiR~2 = 0.956 the mirrors reflectivity, the parameters a r = 150 and a K = 25000 are 
related to the reflectivity of the window W, and are chosen to fit the experimental results. 



In Fig. [| (b) and (d) the numerically generated intensity profiles are compared with their 
experimental counterparts reported in Fig|] (a) and (c), respectively. It can be observed that 
the intensity profiles show antiphase oscillations for both polarization components, while the 
total intensity remains unmodulated. 

When 5=0, the polarization oscillations remain undamped for any degree of anisotropy. 
For 5^0, the oscillations are still undamped only in perfect cavity symmetry conditions 
a = (3 = 0. In Fig.[| (d) it can be seen that for 5 = 0.05, a detuning or losses anisotropy 
of 0.5 % is sufficient to damp the oscillation in a few hundred microseconds as observed 
in the experiment. In the experimental system unavoidable residual anisotropies break the 
cylindrical symmetry, and therefore the polarization oscillations are always damped. 

Choosing the resonant condition 5 = 0, we observe that the frequency of the oscillations 
presents the same dependence on the total losses and pump strength found in the experi- 
ment. In Fig. |6| we compare the experimental and numerical angular dependences of the 
polarization oscillation frequency (Fig. |] (a),(c)) and the total intensity (Fig. | (b),(d)), 
showing good agreements. 

IV. CONCLUSIONS 

In this work we have reported an experimental and theoretical study of the transient 
polarization dynamics of a quasi-isotropic CO2 laser. A competition between two modes 
with orthogonal polarization directions is observed. The competition manifests itself as an 
oscillation on the polarized laser field, whose features depend on the cavity parameters. 

To explain the observed transient polarization alternation, both intrinsic and extrinsic 
anisotropies must be considered; the intrinsic ones provide the alternation and the extrinsic 
ones provide damping and residual ripple of the total intensity. The numerical results fit 
well with the experimental observations, pointing out the importance of matter variables in 
the dynamics of a quasi-isotropic laser. 

The authors are grateful to F.T. Arecchi for fruitful discussions. 
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FIG. 1: Experimental setup: Mi: total reflecting flat mirror, W: additional intracavity window, 
M2: outcoupler mirror, V: Vertical axis of the cavity, H: Horizontal axis of the cavity, Pzt: piezo 
electro translator, Pol: wire grid polarizer, Di: fast detector for the vertical component, D2: fast 
detector for the horizontal component 
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FIG. 2: Experimental intensity of both polarization components along a free spectral range (Black 
triangles vertical polarization direction, gray squares horizontal polarization direction) (a) quasi 
isotropic condition 8 = (b) 8 = 6° (c) 8 = 9°. It is to be noted that a Pzt adjustement is required 
to recover the resonant condition for each 8. 



10 





0.2 0.3 

Time (sec) 



FIG. 3: Intensity of vertical polarization component when the cavity detuning is chosen at the 
transition point near the resonance. 
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FIG. 4: Experimental time intensity profiles for 9 = and slightly different detuning condition : 
(a) resonance, (c) out of resonance. Numerical generated intensity profiles for 0=0, a = (3 = 0.01 
and : (b) 8 = , (d) 5 = 0.05. 



12 




FIG. 5: Numerical intensity profiles for 5 = 0.05, 9=0 and :(a) a = (3=0.005, (b) a = (3=0.0 and 
(c) a = (3=0.005, C=0. 
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FIG. 6: Experimental angular dependences of: (a) polarization oscillation frequency and (b) total 
intensity. Numerical dependences (5 = 0, a=0.01 /3=0.01): (c) polarization oscillation frequency 
and (d) total intensity. 
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